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We present measurements of CP- violating asymmetries in the decay i?" — > of (1260) vr^ with 
af (1260) — > TT^TT^TT*. The data sample corresponds to 384 x 10^ BB pairs collected with the 
BABAR detector at the PEP-II asymmetric S-factory at SLAC. We measure the CP-violating 
asymmetry A'^p = —0.07 ± 0.07 ± 0.02, the mixing-induced CP violation parameter Sai-w = 
0.37 ± 0.21 ± 0.07, the direct CP violation parameter Cai^ = -0.10 ± 0.15 ± 0.09, and the pa- 
rameters ACain = 0.26 ± 0.15 ± 0.07 and ASa^^ = -0.14 ± 0.21 ± 0.06. From these measured 
quantities we determine the angle a^B = 78.6° ± 7.3°. 

PACS numbers: 13.25.Hw, 12.15.Hh, ll.30.Er 



The angle a = arg [-y.^y/^/V^^F^] of the unitar- 
ity triangle of the Cabibbo-Kobayashi-Maskawa (CKM) 
quark-mixing matrix f]\ has recently been measured 
by the BaBar and Belle Collaborations from time- 
dependent CP asymmetries in the decays to tt+tt" 
0, p±7r=F 0, and p+p' The _decay 5° to oitt Q 
proceeds dominantly through the b uud process in 
the same way as the previously studied modes. How- 
ever, due to the presence of additional loop contribu- 
tions, these measurements determine an effective value 
aeff, rather than a itself. This obstacle can be overcome 
using isospin symmetry 0, with bounds to Aa = a — acS 
determined using either an isospin analysis Q or broken 
SU(3) flavor symmetry Q. Because it has the smallest 
contribution from loop diagrams, the B^ P~ decay 
currently allows the most precise single determination of 
a 0. The BaBAR collaboration recently reported the 
observation of B° — > af tt^ where the angle aes can 
be determined by measuring time-dependent CP asym- 
metries 11, l3|. The state afir^, like p^tt^, is not a CP 
eigenstate and four flavor-charge configurations must be 
considered (B"(_B°) afir^). Theoretical bounds on 
Aa in these decay modes based on SU(3) flavor symme- 
try have been derived in Ref. [l^ . 



In this Letter we report measurements of the CP pa- 
rameters in the decay B^ —> af tt^ with af — > tt^tt^tt^. 
The analysis is done in the quasi-two-body approxima- 
tion [13l_. The data were collected with the BaBAR de- 
tector ^] at the PEP-II asymmetric e+e~ collider ^Td\. 
An integrated Imninosity of 349 fb~^, corresponding to 
384 ± 4 million BB pairs, was recorded near the T(45') 
resonance ( "on-resonance" ) at a center-of-mass (CM) en- 
ergy y/s = 10.58 GeV. An additional 37 fb~^ were taken 
about 40 MeV below this energy ( "off- resonance" ) for 
the study of continuum background in which a charm or 
lighter quark pair is produced. 

From a candidate BB pair we reconstruct a B^ de- 
caying into the final state / = aitt (B^^j^). We also 
reconstruct the vertex of the other B meson (-Bj'ag) and 
identify its flavor. The difference At = tain — itag of 
the proper decay times of the reconstructed and tag B 
mesons, respectively, is obtained from the measured dis- 
tance between the B^^j^ and B^^^ decay vertices and from 
the boost (/37 = 0.56) of the e+e~ system. The At dis- 



tributions are given ^12] by: 



"^"^^ ^ l-QtagAiU- 



)tag(l - 2w) 



{Sa^TT ± A^ai^) sin(Am<jAi) - 
(Ca,^ ± ACa.n) cos{AmdAt) 



where Qtag = 1(^1) when the tagging meson B^^^ is a 
B^{B'^)^ t is the mean _B° lifetime, Amd is the mass dif- 
ference between the two B^ mass eigenstates, and the 
mistag parameters w and Aw are the average and dif- 
ference, respectively, of the probabilities that a true 5° 
is incorrectly tagged as a B^ or vice versa. The time- 
and flavor-integrated charge asymmetry A'^p measures 
direct CP violation. The quantities Sai-n and Ca^r: pa- 
rameterize the mixing-induced CP violation related to 
the angle a, and flavor-dependent direct CP violation, re- 
spectively. The parameter ACa^Tr describes the asymme- 
try between the rates T{B^ ^ a'^'^^) + r(B° a^TT+) 
and T{B° aj"7r+) + r{W atn^), while ASa,^ is 
related to the strong phase difference between the am- 
plitudes contributing to B^ aiir decays. The parame- 
ters ACanr and ASa^r are insensitive to CP violation. 
The flavor-tagging algorithm uses six mutually exclu- 
sive categories. Its analyzing power is measured to be 
(30.4 ±0.3)% [B]. 

Charged particles are detected and their momenta 
measured by the combination of a silicon vertex tracker 
(SVT), consisting of five layers of double-sided silicon 
detectors, and a 40-layer central drift chamber, both op- 
erating in the 1.5-T magnetic field of a superconducting 
solenoid. Charged-particle identification (PID) is pro- 
vided by the average energy loss (dE/dx) in the track- 
ing devices and by an internally reflecting ring-imaging 
Cherenkov detector (DIRC) covering the central region. 
Separation between pions and kaons is achieved at the 
level of four standard deviations (a) for momenta below 
3 GeV, decreasing to 2.5 a at 4 GeV. 

Full Monte Carlo (MC) simulations [i3| of the signal 
decay modes, continuum, and BB backgrounds are used 
to establish the event selection criteria. The MC signal 
events are simulated as 5° decays to aiir with ai ~> pir. 
For the ai meson parameters we take the mass mo — 1230 
MeV and the width Fq = 400 MeV [11, Ei. 
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We reconstruct the decay ai —^ Stt with the follow- 
ing requirement on the invariant mass: 0.87 < < 
1.8 GeV. The intermediate dipion state is reconstructed 
with an invariant mass between 0.51 and 1.1 GeV. We 
impose several PID requirements to ensure the iden- 
tity of the signal pions. For the decay pion coming 
from the B meson we require the measured Cherenkov 
angle to be within —2 a and +5 a from the expected 
value for a pion. This requirement removes 98.6% of 
the background from aiK. A B candidate is char- 
acterized kinematically by the energy-substituted mass 
TOES = v^(s/2 -I- po • pb)V-E'o - Pb a-nd energy differ- 
ence AE — Eg — where the subscripts and B 
refer to the initial T(4S') and to the B candidate in the 
laboratory frame, respectively, and the asterisk denotes 
the CM frame. The resolutions in toes and in Ai? are 
about 3.0 MeV and 20 McV respectively. We require 
\AE\ < 0.1 GeV and 5.25 < toes < 5.29 GeV. To reduce 
the number of false _B-meson candidates we require that 
the probability of the B vertex fit be greater than 0.01. 
The absolute value of the cosine of the angle between the 
direction of the tt meson from ai — > pn with respect to 
the flight direction of the B in the ai meson rest frame 
is required to be less than 0.85 to suppress combinatorial 
background. The distribution of this variable is uniform 
for signal and peaks near unity for this background. 

To reject continuum background, we use the angle 9t 
between the thrust axis of the B candidate and that of 
the rest of the tracks and neutral clusters in the event, 
calculated in the center-of-mass frame. The distribu- 
tion of cos 6't is sharply peaked near ±1 for qq can- 
didates, which have a jet-like topology, and is nearly 
uniform for the isotropic _B-meson decays. We require 
I cos 6't I < 0.65. We discriminate further against qq 
background with a Fisher discriminant that combines 
several variables that characterize the production dynam- 
ics and energy flow in the event [20| . The remaining con- 
tinuum background is modeled from off-resonance data. 

We use MC simulations of B'^B^ and B^B^ decays 
to look for BB backgrounds, which can come from 
B decays with or without charmed particles in the fi- 
nal state. Neutral and charged D mesons may con- 
tribute to background through particle mis-identification 
or mis-reconstruction. We remove any combinations of 
the decay products, including possible additional tt*^, 
with invariant mass consistent with nominal mass val- 
ues for K^TT^TT^ or 7r± and K^t:^ or 
K^TT^TT^. The decay mode B° a^(1320)7r=F has the 
same final-state particles as the signal. We suppress this 
decay with the angular variable 7i , defined as the co- 
sine of the angle between the normal to the plane of the 
Stt resonance and the flight direction of the primary pion 
from B meson evaluated in the Stt resonance rest frame. 
Since the ai and a2(1320) mesons have spins of 1 and 2 
respectively, the distributions of the variable Ti, for these 
two resonances differ. We require \H\ < 0.62. 



The average number of candidates found per selected 
event is 1.32. In the case of events with multiple can- 
didates we choose the candidate with the best B- vertex 
fit probability. From simulated signal events we find that 
this algorithm selects the correct candidate in about 92% 
of the events containing multiple candidates, and intro- 
duces negligible bias. 

We obtain the CP parameters and signal yield from 
an unbinned extended maximum likelihood (ML) fit 
with the input observables AE, toes, fnai, 
and At. We have six fit components in the likeli- 
hood: signal, charm and charmless BB background, 
B^ — > a^(1320) TT^, continuum qq background, and non- 
resonant piTTT. The charmless component also includes 
candidates that were incorrectly reconstructed from par- 
ticles in events that contain a true signal candidate. 
Based on measurements of branching fractions for sim- 
ilar charmless decays, we assume 

(2 ± 2) X 10~^, which corresponds to 19 expected events 
in the ML fit sample. This yield is fixed in the fit and a 
systematic uncertainty is assigned to the final results. 

The total probability density function (PDF) for the 
component j and tagging category c in the event i, Vj^^, 
is written as a product of the PDFs of the discriminating 
variables used in the fit. The factored form of the PDF 
is a good approximation since linear correlations among 
observables are below 10%. The systematic uncertainty 
from residual correlations is taken into account in the fit 
bias. We write the extended likelihood function for all 
events as 



3,c 



(2) 



where Uj is the yield of events of component j, fj^c is 
the fraction of events of component j for each category 
c. Tic = fj,cfij is the number of events found by the 
fitter for category c, and Nc is the number of events of 
category c in the sample. We fix c to fBtv^^,c, the val- 
ues measured with a large sample of fully reconstructed 
B^ decays into flavor eigenstates (Bflav sample) [2l[, for 
the signal, /Otttt, and a^(1320) Tr'f' flt components. 

We fix fjc to values obtained with MC events for the 
charmless and charm fit components and allow it to vary 
for the qq component. 

The PDF 7'sig(Ai, tTAt;c), for each category c, is the 
convolution of F{At] c) (Eq. [1]) with the signal resolu- 
tion function (sum of three Gaussians) determined from 
the -Bflav sample. The Ai resolution functions for all the 
other fit components are also modeled with the sum of 
three Gaussians. For charmless, B^ a^(1320)7rT, and 
pTTTT components in the nominal fit to the data we as- 
sume S* = 0, C = 0, AS* = 0, and AC = and we vary 
these parameters when evaluating systematic uncertain- 
ties on final results. We use an effective B lifetime for the 
charmless component as obtained from a fit to MC signal 
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events. The continuum (charm) At distributions are pa- 
rameterized as sums of three Gaussians with parameters 
determined from a fit to off-resonance (MC) events. 

The PDF of the invariant mass of the ai meson in sig- 
nal events is parameterized as a relativistic Breit-Wigner 
hne-shape with a mass-dependent width that takes into 
account the effect of the mass-dependent p width [2^ . 
The PDF of the invariant mass of the 02 (1320) meson 
is parameterized by a sum of three Gaussian function 
distributions. The mss and AE distributions for signal 
are parameterized as a sum of two Gaussian distribu- 
tions. The AE distribution for continuum background 
is parameterized by a hnear function, and the combina- 
torial background in rriEs is described by a phase-space- 
motivated empirical function [i^. We model the Fisher 
distribution J- using a Gaussian function with different 
widths above and below the mean. The A distributions 
are modeled using polynomials. 

The PDF parameters are determined from MC simu- 
lated events with the exception of the continuum back- 
ground, where we use off-resonance data, and of the sig- 
nal resolution function, where we use the Bflav sample. 
Large data control samples of B decays to charmed final 
states of similar topology are used to verify the simulated 
resolutions in toes and AE. Where the control samples 
reveal differences between data and MC in mass and en- 
ergy resolution, we shift or scale the resolution used in 
the likelihood fits. 

We test and calibrate the fitting procedure by apply- 
ing it to ensembles of simulated qq experiments drawn 
from the PDF, into which we have embedded the ex- 
pected number of signal, charmless, i?° — > a^(1320)7rT, 
the charm, and the piriT events randomly extracted from 
the fully simulated MC samples. The measured quan- 
tities S'aiTT, Ca^TT, A^Q^^, ACai,r, and A^p have been 
corrected for the fit biases and a systematic uncertainty 
equal to half of the bias found in MC simulations is as- 
signed on the final results. 

In the fit there are 35 free parameters, including SaiTr, 
CaiTT, ASaiTT, ACaiTT, the charge asymmetries for sig- 
nal and continuum background, five yields, the signal 
fli width, eleven parameters determining the shape of 
the combinatorial background, and 12 tagging efficien- 
cies for the continuum. The main contributions to the 
systematic error on the signal parameters are summa- 
rized in Table |T1 We have studied systematic uncertain- 
ties arising from several sources: variation of the sig- 
nal PDF shape parameters within their errors; model- 
ing of the signal At distribution; tagging efficiency and 
mistag rates determined from the Bflav sample 21|; un- 
certainties in Amd and r [18]; uncertainty in fit bias; 
uncertainty due to CP violation present in the BB back- 
ground, the a^(1320)7r^ CP violation; uncertainty due 



to the interference between 



TT^ and other An 



TABLE I: Summary of the systematic uncertainties (in units 
of 10"^). 





c 

Jain 




A C 

i^Jain 




^aln 


PDF parameterization 


4.8 


5.3 


3.3 


5.3 


1.5 


Signal At model 


0.2 


0.2 


0.3 


0.1 


0.0 


Tagging and mistag 


0.3 


0.2 


0.2 


0.4 


0.1 


Anid and r 


0.0 


0.2 


0.3 


0.1 


0.0 


Fit bias 


0.8 


0.2 


0.8 


1.0 


0.3 


BB CP violation 


4.1 


4.3 


4.2 


4.0 


0.5 


aJ(1320)7r'F + interf. 


2.8 


4.5 


3.2 


0.6 


0.2 


DCS decays 


0.8 


2.2 


0.0 


2.2 


0.1 


SVT alignment 


1.0 


0.6 


1.0 


0.6 


0.0 


Particle ID 


0.1 


0.1 


0.1 


0.1 


0.2 


Total 


7.0 


8.5 


6.4 


7.1 


1.6 



for some tag-side B decays 24]; SVT alignment; and 
the particle identification algorithm. We allow for a CP 
asymmetry up to 20% in B decays to charmless final 
states, and up to 50% in B decays to a2(1320)7r. 

From the fit to a sample of 29300 events, we obtain 
a signal yield of 608 ± 53, of which 461 ± 46 have their 
fiavor identified and are used to measure the following ad- 
ditional parameters: Sa^n ~ 0.37 ± 0.21 ± 0.07, ASa^Tv — 
-0.14±0.21±0.06, Ca,n = -0.10±0.15±0.09, ACa^^ = 
0.26 ± 0.15 ± 0.07, Ac'p = -0.07 ± 0.07 ± 0.02. Linear 
correlations between these fit parameters are small. 

The angle can be defined 12] as: 



arcsm 



Sa 1 7r 



^l-iCa,.+ACa,.y, 



(3) 



arcsm 



-\/l — {CaiiT — ACaiir)^ 



final states have been estimated with MC simulations; 
doubly-Cabibbo-suppressed (DCS) b — > ucd amplitude 



Using the measured CP parameters, we determine the an- 
gle fteff and one of the four solutions, tteff = 78.6° ± 7.3°, 
is compatible with the result of SM-based fits. Using the 
published branching fraction [loj and adding statistical 
and systematic errors in quadrature, we obtain also the 
following values for the flavor-charge branching fractions 
^ (in units of 10-^): B{B° a+T^) = 17.9 ± 4.8, 
S(SO ^ a-n+)^ 11.4 ± 4.7, B{B° ^ a+Tr") = 
13.0 ± 4.3, and B{B° a^TT+) = 24.2 ± 5.8. 

Figure [1] shows distributions of mEs and AE, en- 
hanced in signal content by requirements on the signal-to- 
continuum likelihood ratios using all discriminating vari- 
ables other than the one plotted. Figure [2] gives the At 
projections and asymmetry for flavor tagged events. 

In summary, we have measured the CP- violating asym- 
metries in i?" af (1260) tt^ decays and determined the 
angle acs- We do not find evidence for direct or mixing- 
induced CP violation in these decays. Once measure- 
ments of branching fractions for SU(3)-related decays 
become available, quantitative bounds on Aa obtained 
with the method of Ref. will provide significant con- 
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FIG. 1: Projections of a) AE, b) rriES. Points represent on- 
resonance data, dotted lines the sum of all backgrounds, and 
solid lines the full fit function. These plots are made with a 
cut on the signal likelihood. 
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FIG. 2: Projections onto At of the data (points) for a) and 
b) -B" tags, showing the fit function (solid line), and the back- 
ground function (dotted line), and c) the asymmetry between 
B° and B° tags. 



straints on the angle a through the measurement of a^s 
in B° a]^(1260)7rT decays. 
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